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abstract 


An  experimental  study  of  the  forced  convection  heat  loss 
from  long  thin  cylinders  placed  normal  to  an  Incompressible  flow 
has  been  carried  out.  Tests  were  made  in  three  different  gases 
(air,  argon,  helium)  at  atmospheric  pressure  and  In  air  over  a 
pressure  range  0.2  to  8  atm.  Particular  attention  was  given  to 
determining  the  effect  of  cylinder  temperature  on  the  heat  loss. 

It  was  found  that  the  "constants"  employed  in  the  usual  empirical 
representation  of  this  heat  loss  vary  with  cylinder  temperature  in 
a  way  dependent  on  Knudsen  number,  results  in  helium  giving  a 
marked  indication  of  this.  A  velocity  fluctuation  sensitivity 
equation  was  derived  which  includes  the  effects  of  this  temperature 

variation. 
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HOT-WIRE  HEAT  LOSS  AND  FLUCTUATION  SENSITIVITY 
FOR  INCOMPRESSIBLE  FLOW 


I.  INTRODUCTION 

Although  the  hot-wire  anemometer  has  been  employed  both  as  an 
AC  and  a  DC  instrument  tor  more  than  thirty  years,  the  general  relation 
between  the  heat  loss  from  the  wire  and  the  properties  of  the  stream 
passing  over  it  remains  uncertain.  Neither  precise  theoretical  formulae 
nor  general  experimental  data  for  the  total  heat  loss  from  long  cylinders 
placed  normal  to  a  stream  at  low  but  finite  Reynolds  numbers  are 
available.  The  few  theoretical  treatments  of  the  problem  suffer  from 
the  imposition  of  unrealistic  simplifications;  the  experimental  studies 
from  either  lack  of  precision  or  insufficient  generality.  Only  one 
valid  limiting  case  has  been  solved  theoretically.  It  was  given  In 
simplest  form  comparatively  recently  by  Cole  and  Roshko  (1954).  In 
this  solution,  Oseen  flow  was  assumed;  the  momentum  and  energy  equ¬ 
ations  were  thereby  rendered  independent  and  the  latter  was  solved 
for  the  total  heat  loss  from  an  infinite  cylinder.*  The  theory  is 
valid  only  for  the  limiting  case  R  -•  0;  in  addition,  since  a  constant 
thermal  conduct i vh ty  was  assumed,  it  applies  only  for  small,  temperature 
differences  between  the  cylinder  and  the  stream.  Experiment  (Col  I  is  and 
Williams,  1959)  indicates  an  approach  to  the  theoretical  result  for 
very  small  Reynolds  numbers.  An  extension  of  this  theory  to  include 
temperature-jump  effects  for  near-cont 1 nuum  flow  about  the  cylinder 
has  been  given  (Levey,  1959)  but  its  verification  requires  a  knowledge 


*  For  reference  to  other  Oseen-type  solutions,  see  Corrsin  (1961). 


of  the  thermal  accomodation  coefficient  which  Is  generally 
difficult  to  obtain. 

For  Reynolds  numbers  of  the  order  of  unity  or  greater, 
the  form  of  the  heat  loss  law  has  conventionally  been  assumed  to  be 
that  of  the  empirical  King's  law  (King,  1914),  which,  expressed  in 
non-dimensional  form,  is 

Nu  =  A  +  B  R 

where  A  and  B  are  considered  as  constants,  TNu  is  the  wire  Nusselt 
number,  R  the  Reynolds  number.  A  variety  of  reference  temperatures 
upon  which  to  base  the  conductivity,  density,  and  viscosity  have  been 
employed  in  this  law.  While  the  values  of  the  quantities  A  and  B  re¬ 
ported  by  various  investigators  are  diverse,  usage  of  the  square  root 
of  the  Reynolds  number  is  fairly  uniform  (Morkovin,  1956,  p.  22; 
Kovasznay,  1954,  p.  241)*.  A  notable  exception  is  the  careful  series 
of  tests  reported  recently  by  Col  I  is  and  Williams  (1959).  There,  a 
Reynolds  number  exponent  of  0.45  was  found  to  fit  the  data  best  in 
the  range  0.02  <  R  <  44. 

Although  there  has  been  considerable  experimental  emphasis 
recently  on  hot-wire  heat  loss  in  compressible  flow,  the  many  details 
of  the  incompressible  case  which  have  not  yet  been  adequately  explored 
leave  ample  room  for  study.  The  effect  of  wire  temperature  on  the 
heat  loss  characteristics  is  uncertain,  heat  loss  measurements  in 
gases  other  than  air  are  scarce,  and  the  effect  of  Knudsen  number  on 
incompressible  heat  loss  is  virtually  unexplored.  A  knowledge  of 

*  The  explanation  often  given  for  this  usage  is  the  square-root 
dependence  with  Reynolds  number  of  the  heat  loss  in  boundary  layer 
problems;  such  a  rationalization  is,  however,  difficult  to  justify 
for  the  low  Reynolds  numbers  commonly  employed  in  hot-wire 
anennometry . 
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these  effects  is  not  only  of  interest  in  itself  but  is  crucial  for 
precise  operation  of  the  hot-wire  as  a  fluctuation-sensing  device. 

In  this  paper,  a  series  of  heat  loss  measurements  conducted  in  three 
gases  (air,  argon,  helium)  and  at  several  ambient  density  levels, 
with  special  attention  to  the  effect  of  wire  temperature,  is  first 
presented.  An  analysis  of  the  fluctuation  sensitivity  of  the  hot¬ 
wire  based  on  these  heat  loss  results  follows. 

In  view  of  the  divergence  in  opinion  over  the  proper  Reynolds 
number  exponent  to  use,  the  more  general  heat  loss  law, 

I )  Nu  =  A  +  BR" 

has  been  assumed  for  fitting  the  data  obtained.  The  wires  were  oper¬ 
ated  at  a  series  of  wire  temperatures  and  the  possible  temperature 
dependence  of  n  as  we  I  I  as  A  and  B  was  explored  by  least— squares 
fitting  to  the  general  form  (I)  of  the  set  of  Nu  vs.  R  data  obtained 
at  each  wire  temperature.  The  reference  temperature  upon  which  to 
base  the  fluid  properties  was  arbitrarily  taken  to  be  the  usual  film 
temperature  (mean  of  wire  and  free  stream  temperatures)  since  there 
is  some  evidence  that  the  temperature  effect  is  reduced  by  using  this 
base.  Conversion  to  other  reference  temperatures  is,  however,  a 
simple  matter. 
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I  I .  APPARATUS 
A.  Tunnel 

For  all  the  measurements  reported  herein,  the  variable  density 
low  speed  blow  down  tunnel  of  the  Gas  Dynamics  Laboratory  was  employed. 
This  tunnel  operates  over  a  range  of  0. 1  to  10  atmospheres  ambient 
pressure  level.  The  piping  arrangement  for  the  supply  of  air  to  the 
tunnel  allows  use  of  either  banks  of  3,000  psi  air  stored  external  to 
the  laboratory  or  of  an  auxiliary  regulated  air  supply  at  a  consider¬ 
ably  lower  pressure.  For  operation  with  gases  other  than  air,  a  simple 
T-connection  permits  the  air  supply  to  be  shut  off  and  other  gases  to 
be  introduced.  These  gases  were  stored  in  external  pressure  regulated 
trailers;  they  had  a  specified  purity  of  99.995?  and  a  moisture  content 
of  less  than  0.02  mg/ 1  iter.  No  information  is  available  on  the  proper¬ 
ties  of  the  air  used;  however,  some  drying  of  the  air  is  accomplished 
prior  to  its  storage  in  the  external  tanks.  In  the  case  of  the  low 
density  runs,  the  tunnel  exit  was  connected  to  the  air  ejector  system 
of  the  laboratory. 

The  free  stream  turbulence  level  of  the  tunnel  was  suffi¬ 
ciently  low  so  as  not  to  influence  any  of  the  heat  loss  results.  A 
porous  steel  filter  of  35  micron  pore  size  was  installed  at  the  tunnel 
inlet;  it  virtually  eliminated  the  problems  of  hot  wire  breakage  and 
dirt  accumulation.  To  measure  high  static  pressure,  a  Bourdon-type 
absolute  gage  with  range  0-200  inches  mercury  and  scale  divisions  of 
0.1  inches  mercury  was  used.  This  gage  had  been  previously  calibrated 
and  the  appropriate  corrections  were  applied  in  each  case  to  the 

For  still  higher  pressures,  a  gage  with  range 


indicated  pressures. 


0-200  psi  and  scale  divisions  of  0.5  psi  was  employed.  In  the  low 
pressure  tests,  a  bank  of  U-tube  mercury  manometers  was  used.  One 
end  of  the  manometer  tubes  was  left  open  to  the  atmosphere,  and  the 
atmospheric  pressure  recorded  from  a  wall-mounted  mercury  barometer. 
Since  the  running  time  was  limited  in  these  tests  due  to  the  high  mass 
flow  rate  of  the  ejector,  all  the  pressure  data  were  taken  on  film  by 
means  of  a  recording  camera.  A  total  head  pitot-probe  in  conjunction 
with  an  inclined-tube  manometer  was  employed  to  measure  the  mean  flow 
velocity  in  the  tunnel.  The  manometer  had  a  range  of  0-2  inches  of 
manometer  fluid  with  0.01-inch  scale  divisions. 

B.  Hot-Wires 

The  physical  and  geometric  properties  of  the  wires  employed 

are  summarized  in  the  following  table: 

Table  I:  Wire  Properties 


D i ameter , 
Wire  d 

Material  (in.xlO  ) 

Fi  rst 
coeff.  of 
Resistance 

(°C“^I0^) 

Second 
coeff.  of 
Resistance 
0o 

(°C"2x 10^) 

Conduc¬ 

tivity 

k 

(watts/ cm°C) 

Heat 

Re lat i ve  Capac i ty 
Density  Cw 

Sw  (cal/gm) 

Length 

Jb 

(in.) 

90Pt-10Rh 

2.5 

1 .7 

-l.0(?) 

0.30 

21  .4 

0.035 

0.08 

Pt  1.0 

3.6 

-7.0(?) 

0.75 

21  .4 

0.035 

0.08 

The  data  in  this  table  are  a  combination  of  nominal  values  given 
by  the  manufacturer  (composition,  diameter),  reference  values  taken  from 
either  handbooks  or  other  sources  (Po^  cw) ,  and  measured  values 

averaged  over  al  I  the  wires  used  (£,  C3q)  . 

The  manufacturer  of  a  1 1  the  wires  was  the  Sigmund  Cohn  Company, 
Mount  Vernon,  New  York.  The  nominal  value  of  diameter  given  by  the 
supplier  has  been  found  by  other  workers  to  be  fairly  accurate. 
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Lawrence  and  Landes  (1952)  report  a  series  of  e lectron-microphotographs 
taken  of  various  samples  supplied  by  the  same  manufacturer.  The  agree¬ 
ment  in  diameter  is  generally  about  5%  with  some  variation  from  sample 
to  sample. 

The  reference  quantities  given  in  the  table  above  have  been 
selected  by  averaging  these  properties  over  the  slightly  conflicting 
values  which  are  reported  from  various  sources.  In  general,  the 
agreement  between  these  sources  is  of  the  order  of  5^.  The  exception 
is  the  value  of  the  "second  temperature  coefficient  of  resistance", 

Pq,  which  occurs  as  the  quadratic  coefficient  in  the  empirical  fitting 
of  wire  resistance  vs.  temperature  given  in  (24).  Since  for  moderate 
temperatures  and  pure  substances,  the  resistance  is  closely  linear  with 
temperature,  the  value  of  is  usually  small  and  is  therefore  diffi¬ 
cult  to  measure  precisely.  In  addition,  Po  seems  to  vary  somewhat  with 
the  temperature-history  of  the  wire  and  with  the  actual  size  of  the 
wire. 

The  values  given  as  measured  in  the  table  are  average  values 
for  the  wires  for  which  such  measurements  were  taken.  The  measurement 
of  Oq  was  accomplished  by  inserting  the  wire  in  an  Insulated  container 
and  recording  the  resistance  as  a  function  of  ambient  temperature. 

The  values  of  I  were  measured  by  a  travelling  microscope. 

The  wire  supports  were  common  steel  sewing  needles.  A  high 
resistance  sealer,  de  Khot i nsky  cement,  was  used  to  fix  the  probe 
leads  and  needles  to  a  ceramic  cylinder.  A  typical  probe  head  is 
shown  in  Figure  I.  Two  of  the  actual  wires  employed  are  shown  in 
Figure  2.  The  wires  were  etched  from  larger  Wollaston  wires  which 
had  been  soldered  to  the  needle  tips. 


6 


111.  EXPERIMENTAL  PROCEDURE 
A.  Cold  Resistance 

As  is  customary,  the  wire  comprised  one  arm  of  a  Wheatstone 
bridge.  The  heated  wire  resistance  was  set  to  a  predetermined  value 
by  adjusting  the  bridge  current  and  the  current  was  measured  by 
applying  a  potentiometer  across  a  precision  resistor  in  series  with 
the  wire. 

In  order  to  set  the  wire  temperature  at  a  fixed  value,  it  is 
first  necessary  to  know  the  resistance  of  the  unheated  wire  at  the 
ambient  temperature  of  the  tunnel.  Since  the  resistance  is  measured 
by  a  bridge  arrangement,  some  current  must  be  passed  through  the  wire. 
This  current  will,  of  course,  induce  some  amount  of  heating  of  the  wire 
and  consequently  a  method  must  be  devised  to  extrapolate  the  resistance 
measured  back  to  a  condition  of  zero  heating.  The  conventional  method 
used  i s  as  fo I  lows : 

1)  With  a  flow  of  fixed  velocity  over  the  wire,  a  small 
known  current,  i^  (say  2  ma) ,  is  passed  through  the  wire, 
and  the  wire  resistance,  R|,  is  measured. 

2)  A  current  I2  =  V  2  i|  's  then  passed  through  the 
wire  and  the  wire  resistance  R2  measured. 

3)  The  unheated  wire  resistance  may  then  be  computed  as 

2)  Rq  =  2  RfRz 

The  relation  (2)  may  be  easily  shown.  The  heat  loss  law  (I) 
for  wire  temperatures  close  to  ambient  is,  in  dimensional  form, 

~  (R-Rq)  CI+cU") 
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where  c  is  a  constant  and  U  Is  the  stream  velocity.  Thus  for 
two  different  currents  at  fixed  velocity, 

|2|  Rj-Rq 

1^2  R2-R0 

If  I2  =V2  l|j  then  (2)  Is  immediately  obtained. 

It  is,  of  course,  necessary  to  subtract  the  lead  resistance 
from  the  measured  values.  The  lead  resistance  was  obtained  by  inserting 
the  unetched  hot-wire  probe  into  the  circuit  and  taking  the  resistance 
measured  to  be  the  lead  resistance.  This  procedure  was  accurate  since 
the  contribution  of  the  resistance  of  the  unetched  portion  of  the  wire 

to  the  total  lead  resistance  was  negligible. 

In  many  cases,  the  hot-wire  filaments  underwent  considerable 
cold  resistance  changes  during  initial  testing.  This  phenomenon  has 
been  noted  by  many  investigators.  The  initial  rapid  decrement  observed 
in  the  resistance  is  probably  due  to  the  annealing  effect  of  high  wire 
temperatures  which  relieves  the  cold-working  Imposed  by  the  drawing 
process  used  In  manufacture.  Subsequently,  however,  the  wire 
resistance  usually  slowly  increased.  A  satisfactory  explanation  of 
this  increase  has  not  been  given.  it  may  be  a  result  of  the  purging 
of  adsorbed  gas  or  other  impurities  in  the  wire.  To  stabilize  the 
cold  resistance,  the  wire  was  operated  at  the  maximum  temperature  to 
be  used  in  the  measurements  and  the  heated  resistance  monitored. 

After  some  length  of  time,  which  varied  from  minutes  to  hours  among 
wire  samples,  the  resistance  usually  stabilized  and  in  subsequent 
tests  remained  nearly  constant. 
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B.  Heat  Loss 


To  measure  the  flow  velocity,  a  pitot  tube  was  stationed 
close  to  the  wire.  At  a  constant  wire  temperature,  the  mass  flow 
over  the  wire  was  varied  and  the  wire  current  and  the  dynamic  pressure 
measured.  Another  static  tap  connected  to  an  appropriate  manometer 
yielded  the  absolute  pressure  level.  In  the  low  density  tests,  the 
pressure  level  of  the  tunnel  was  held  fixed  and  the  velocity  alone 
varied.  This  was  made  possible  by  means  of  a  variable-area  choking 
nozzle  at  the  tunnel  exit.  These  tests  were  therefore  performed  at 
constant  wire  Knudsen  number. 

All  of  the  tests  were  conducted  over  a  range  of  wire  Reynolds 
numbers  (where  the  viscosity  g,  ,  and  the  density  p,  are  based  on  film 
temperature)  of  about 

0.4  <  R  <  4 

This  restricted  range  was  used  to  avoid  possible  changes  in  the  heat 
loss  characteristics  due  to  the  appearance  of  discrete  vortices  on  the 
leeward  side  of  the  wire.  The  Reynolds  number  at  which  these  form  is 
general ly  bel ieved  to  be  about  4. 

Considerable  care  was  exercised  to  insure  that  the  ambient 
stream  temperature  remained  constant  during  the  heat  loss  tests. 
Variations  of  as  much  as  l°F.  in  this  temperature  during  a  run  were 
uncommon.  All  of  the  electrical  and  aerodynamic  quantities  measured 
were  checked  from  time  to  time  by  independent  measurement  using  alternate 
equipment;  excellent  agreement  was  always  found. 
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IV..  DATA  REDUCT  I  ON 
A.  General  Remarks 

Preliminary  evaluation  of  the  data  indicated  that  at  least 
8  evenly  distributed  data  points  for  each  test  were  required,  and 
preferably  10  or  12,  for  consistency.  Consequently,  all  tests  having 
fewer  points  than  8  were  not  considered  in  the  evaluation  of  the  heat 
loss  data  and  if  any  discrepancies  occurred  among  tests  under  the  same 
operating  conditions,  preference  was  given  to  tests  containing  the  most 
number  of  points.  Discrepancies  among  tests  with  |0  or  more  points 
were  rare  and  usually  traceable  to  some  gross  procedural  error. 

The  reduction  of  all  data  was  carried  out  on  an  IBM  650 
electronic  digital  computer.  The  use  of  such  a  computer  is  virtually 
mandatory  for  least-squares  data  fitting  of  the  numerous  heat  loss 
tests  carried  out. 

For  the  density,  the  perfect  gas  equation  of  state  was  assumed. 
The  mean  free  path  was  computed  from  (Kennard,  1938) 

4)  \  =  U  V  TT/2RgT  , 

(where  u  =  p/p  and  Rg  is  the  gas  constant),  and  was  evaluated  at  the 
wire  temperature  and  in  some  cases  at  the  reference  temperature. 

The  physical  properties  of  the  wires  used  are  given  in  Table  I. 
Since  the  principal  emphasis  In  this  work  is  on  determining  the  temp¬ 
erature  dependence  of  A,  B,  n  as  well  as  on  testing  in  several  gases 
at  various  Knudsen  numbers,  no  special  effort  was  made  toward  fixing 
precisely  the  absolute  values  of  A,  B.  Accuracy  in  these  quantities 
depends  principally  upon  accurate  (and  frequent)  determination  of 
OqRo/j^  wire  diameter.  The  consistency  of  the  results  using 
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averaged  values  for  these  quantities  In  the  data  reduction  is  the 
only  available  method  of  estimating  errors  due  to  their  variation. 

Sample  measurements  of  Oq  agree  well  with  values  reported  elsewhere, 
however,  the  range  of  temperatures  employed  was  insufficient  to  detei 
mine  precisely  the  value  of  the  second  coefficient,  Since  Pq  's" 

at  any  rate  expected  to  be  small  for  Pt-Rh  wires,  it  was  taken  to  be 
zero  in  the  computations.  A  discussion  of  the  possible  error  introduced 
by  this  assumption  is  given  later. 

The  general  definition  of  Nusselt  number  is 


where  the  film  coefficient,  h,  is 


B.  Computation  Formulae 
Nusselt  number: 


5) 


Nu  = 


hd 


h 


q 

Aw  ("Tw'To) 


and  q  is  the  totai  heat  transfer.  Aw  the  wire  surface  area.  Thus 
for  anohmicaily  heated  infinite  cylinder. 


|2rw 

^  Nu  - 

If  the  conductivity  is  evaluated  at  the  film  temperature 

Tw+To 

^ 


then  assuming  a  linear  resistance-temperature  relation. 


7) 


Nu 


|2  l+Bw  OoRo/-^ 


where 


i  +h3^ 


Rw“Ro 


is  the  so-called  "overheating  ratio"  and 


(dk/dT)m 

h  =  — 

2koQo 

For  purposes  of  data  reduction,  it  is  more  convenient  to  use  aw  as  a 
measure  of  wire  temperature  rather  than  the  temperatur'e  loading 
quant i ty 

^w'Tq 

T  =  — -  . 


Corrections:  To  account  for  loss  of  heat  to  the  wire  supports 
by  conduction  the  approximate  relation  for  an  "infinite  wire"  Nusselt 
number , 

.  Nu 

8)  =  -  ■ 

1+1 .  IS/  V  l+ow 

given  by  Kovasznay  (1954)  was  used,  where 

d  /  l^w 

Betchov  (1952)  has  studied  the  end  loss  correction  when  the  quantity 
A  =  A(t)  and  it  may  be  shown  from  his  results  that,  to  order  S,  the 
effect  of  end  loss  and  variable  A  are  independent.  The  correction  to 


12 


order  S  is  given  by  (8)  except  that  the  factor  l.l  Is  replaced  by 
unity.  This  factor  was  used  by  Kovasznay  to  give  a  slightly  better 
approximation  to  the  end  loss  correction  while  retaining  only  the 
first  order  term  In  S. 

A  method  of  computing  the  effect  of  "temperature-jump"  in  near 
continuum  flow^  suggested  by  Qollis  and  Williams  (1959)  and  found  theo¬ 
retically  by  Levey  (1959),  simply  involves  correcting  the  temperature  of 
the  gas  at  the  wire  surface  by  the  amount  predicted  for  the  temperature 
Jump  from  elementary  kinetic  theory.  It  yields  for  the  continuum 
Nusselt  number 


10) 


Nu 

I  “CKf^Nu 


where 

II) 


2(2-f)  Y  I 
f  Y+r  Pr 


and  f,  y,  Pr  are  the  thermal  accomodation  coefficient,  specific  heat 
ratio,  and  Prandt I  number,  respectively.  The  applicability  of  this 
simple  correction  outside  the  range  R  -•  0  remains  to  be  demonstrated; 
however,  for  want  of  a  more  precise  method,  it  was  used  in  the  data 
reduction  here.  Even  assuming  its  validity,  there  is  considerable 
question  about  what  value  to  use  for  the  accomodation  coefficient,  f. 

As  pointed  out  by  Co i I  is  &  Wi I  I iams,  substitution  of  a  typ i cal  va lue 
of  0.9  reported  in  the  literature  for  air  on  platinum  gives  a  value 
of  C  for  air  very  close  to  2  which  was  consequently  employed  in  their 
data  reduction.  The  same  value  was  used  here  for  both  air  and  argon 
on  both  Pt  and  90Pt-l0Rh  wires.  The  case  of  helium  is  discussed  later. 
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Combining  the  two  corrections  above  gives 


Nu 

12)  Nu*  = - -  ,  — 

(l-CKnNu)  ( l+l  •  IS/  V  H«w 

as  the  continuum  infinite-wire  Nusselt  number. 

Reynolds  number:  The  wire  Reynolds  number  is 


R 


pUd 

P> 


Based  on  film  temperature,  this  becomes 

PoUd 

13)  R  = -  (i+h|aw)  (l+b2aw)> 

M-o 


where 

{dii/dT)m  ^  I 

h  I  =  - ,  _  f  h2  =  2a  T 

Knudsen  number;  The  wire  Knudsen  number  is 
X 

K„._ 

and  using  equation  (4), 

Kn  =  V  TT/2Rgf 

Since  non-continuum  effects  probably  depend  on  the  properties  at  the 
wire  rather  than  in  the  external  field,  these  properties  were  usually 
evaluated  at  the  wire  temperature,  giving 

14)  '^n  =  ^  V  TT/2RgTo  ( 1 +2h  ,a„)  (  1 +2h2a«)^ 
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Estimates  of  radiation  and  free  convection  effects  indicated  that 
these  were  entirely  negligible  over  the  range  of  this  experiment. 
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V,  DISCUSSION  OF  HEAT  LOSS  RESULTS 
A.  Air  and  Argon  Results  at  Atmospheric  Pressure 

Sample  data  for  a  IxlO'^  In.  diameter  90Pt-l0Rh  wire  are  shown 
In  Figure  3.  This  heat  loss  measurement  was  carried  out  in  air  at 
atmospheric  pressure  with  a  wire  temperature  of  about  88CPF.  In  the 
upper  diagram,  the  data  are  plotted  In  the  usual  coordinates  used  for 
hot-wire  calibrations.  The  convex  curvature  when  plotted  with  the 
square  root  of  velocity  is  apparent.  From  the  raw  data  of  this  figure, 
the  Nusselt  number  vs.  Reynolds  number  plot  (based  on  film  temperature) 
of  the  lower  figure  was  constructed.  Least-squares  fitting  of  the  data 
Indicated  that  the  proper  Reynolds  number  exponent  for  these  data  was 
0.34.  No  end- loss  or  temperature  jump  corrections  have  been  applied 
to  the  Nusselt  number. 

To  determine  the  variation  of  A,  B,  n  with  wire  temperature, 
a  series  of  heat  loss  measurements  in  air  and  argon  were  carried  out 
with  a  given  wire  over  a  range  of  wire  temperatures.  In  each  case, 
both  the  corrected  and  uncorrected  Nusselt  numbers  were  computed  and 
their  variation  with  Reynolds  number  subjected  to  least-squares  fitting 
The  wire  used  was  90Pt-l0Rh  with  a  diameter  of  IxlO"^  Inches.  The 
variation  of  the  quantities  A  and  B  as  we  I  I  as  the  exponent  n  with 
T  +  I  are  shown  in  Figure  4  for  the  case  wherein  corrections  for 
end- loss  and  temperature  jump  have  been  applied.  The  Knudsen  number 
for  the  latter  correction  was  based  on  wire  temperature.  Included  In 
these  plots  are  also  the  results  from  three  heat  loss  measurements 
carried  out  with  a  different  wire  (designated  wire  II)  of  the  same 
material  and  diameter.  The  solid  lines  represent  the  least-squares 
fitting  of  the  logarithm  of  each  quantity  vs.  Xn(T+l).  Thus  the 
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values  at  the  abscissa  origin  are  the  heat  loss  "constants"  in  the  iimit 
T  -  0.  To  the  right  of  each  curve  are  given  the  logarithmic  slopes 

determined  by  the  fitting  process. 

It  may  be  seen  that  the  variation  of  A  with  wire  temperature  Is 
very  small  and  is  probably  within  the  experimental  accuracy  rvolved  in 
determining  its  variation.  Further,  the  agreement  between  the  actual 
numerical  value  of  A  for  wires  30  and  II  is  good.  This  is  an  indication 
of  the  closeness  in  OoRo/X  between  the  two  wires.  There  is,  however,  a 
significant  difference  in  A  for  air  and  argon. 

The  variation  of  B  with  temperature  has  a  definite  non-zero 
value;  it  appears  to  vary  with  T+i  to  about  the  i/5  power  for  argon  and 
to  the  i/7  power  for  air.  Also  the  values  of  B  are  somewhat  higher  in 
argon  than  in  air.  While  the  dependence  of  B  on  temperature  appears  to 
be  about  the  same  for  wire  il  as  for  wire  30,  the  values  of  B  are 
slightly  higher  in  the  former  case.  This  difference  probably  results 
from  a  difference  in  actual  diameter  of  the  two  wires  (the  nominal  value 
was  used  in  the  computations). 

With  regard  to  the  variation  of  n  with  temperature,  the  vari¬ 
ation  is  small  but  probably  significant.  in  argon,  the  logarithmic 
’slope  is  about  i/iO  and  in  air  about  i/8.  The  numerical  values  of  n 
agree  well  between  wire  30  and  wire  I  I.  in  argon,  however,  an 
exponent  consistently  smaller  than  that  in  air  was  found. 

B.  Results  for  Small  Knudsen  Number  in  Air 

Turning  to  Figure  5,  the  variation  of  the  quantities  A,  B,  n 
were  here  determined  from  data  supplied  by  heat  loss  measurements  at 
high  ambient  density  levels.  in  general,  the  pressure  for  each  of 
these  measurements  varied  from  about  2  atmospheres  to  8  atmospheres. 
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Thus,  although  the  Knudsen  number  based  on  wire  diameter  and  mean 
free  path  of  the  free  stream  varied  during  a  given  test,  in  each  case 
the  pressure  level  was  sufficiently  high  to  render  the  Knudsen  number 
considerably  less  than  0.02  and  hence  the  temperature  Jump  correction 
essentially  negligible  for  the  lx  10"'^  in.  diameter  wire.  Also 
plotted  in  Figure  3  are  values  of  A,  B,  n  obtained  for  a  wire  of 
2.5x10“'^  in.  diameter  at  atmospheric  pressure.  Here  again  the  Jump 
correction  is  small.  Nevertheless,  in  each  case  temperature  jump 
corrections  using  a  Knudsen  number  based  on  wire  temperature  as  well 
as  end  loss  corrections  have  been  applied. 

The  variation  of  A  with  temperature  appears  to  be  somewhat 
larger  than  in  the  previous  case;  it  appears  here  to  vary  with  about 
the  l/4  power  of  wire  temperature  for  the  small  wire.  For  the  large 
wire,  a  similar  variation  with  temperature  appears,  however,  there 
is  a  striking  difference  in  the  magnitude  of  A  between  the  small 
and  large  wires.  No  satisfactory  explanation  for  this  difference 
is  apparent;  it  seems  to  be  well  above  possible  errors  due  to 
incorrect  determination  of  a^R^/ 1. 

The  quantity  B  varies  strongly  with  temperature,  considei — 
ably  more  so  than  in  the  case  at  atmospheric  pressure.  The  variation 
here  with  about  the  l/3  power  of  temperature  is  greater  than  for  any, 
other  quantity.  It  may  be  noticed  that  a  similar  variation  seems  to 
occur  for  the  large  wire  as  well  and  further  that  the  numerical 
values  of  B  agree  well  between  the  two  wires.  The  latter  is  a  good 
indication  that  the  actual  wire  diameters  were  close  to  the  nominal 
ones  used  in  the  computation. 

A  variation  of  n  with  t+I  very  nearly  identical  to  the 
variation  at  atmospheric  pressure  may  be  seen  in  the  lowest  curves 
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of  Figure  5.  The  agreement  between  the  large  wire  and  the  small 
wire  at  high  density  is  also  quite  close. 

C.  Large  Wire  at  Low  Density 

The  final  case  studied  in  air  was  the  heat  loss  from  the 
large  wire  at  low  density.  Here,  each  heat  loss  measurement  was 
carried  out  at  a  fixed  free  stream  density  level,  lower  than  atmos¬ 
pheric,  and  over  a  range  of  wire  temperatures.  It  should  be  pointed 
out  that,  because  of  the  limited  running  time  for  each  heat  loss 
measurement  at  low  density  a  considerable  diminuition  of  experimental 
accuracy  was  inevitable  and  consequently  a  larger  scatter  in  the  data 
I s  to  be  expected.*  Nevertheless,  a  fair  amount  of  consistency  with 
the  previously  discussed  results  was  obtained.  The  variation  of  A, 

B,  n  with  wire  temperature  at  a  free  stream  density  level  of  about 
0.2  atmospheres  is  shown  in  Figure  6.  Values  corrected  for  end  loss 
and  temperature  jump  as  well  as  uncorrected  values  are  plotted 
together.  In  addition,  the  quantities  from  Figure  5  for  the  large 
wire  at  atmospheric  pressure  are  also  plotted.  A  set  of  heat  loss 
measurements  at  about  0.5  atmospheres  with  the  large  wire  was 
carried  out  but  is  not  shown;  the  results  will  be  given  later. 

Considering  first  the  variation  of  A  with  temperature  in 
Figure  6,  it  is  seen  that  the  corrected  data  at  low  pressure  show 
a  fair  agreement  with  the  atmospheric  pressure  data.  To  one  signi¬ 
ficant  figure  (which  in  view  of  the  scatter  is  certainly  the  most 
that  may  be  considered),  the  slope  A  agrees  with  the  value  from 
Figure  5.  That  the  corrections  for  end  loss  and,  particularly,  for 

*  The  weak  dependence  of  A,  B,  n  on  wire  temperature  calls  for 
extremely  accurate  and  consistent  measurements  in  order  to  be 
able  to  detect  this  dependence;  it  is  believed  that  this  fact 
is  principally  responsible  for  the  divergence  of  opinion  in  the 
I iterature  over  the  proper  value  for  these  quantities, 
part  icu lar ly  for  n . 
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i'GmpGrai'ur©  jump  havG  a  s  i  9^  *  ©ttGci*  on  i'hG  slop©  is  ©vidGni* 

from  a  comparison  of  th©  corr©cf©d  and  uncorr©cted  valuGs  of  A. 

Th©  fact  that,  wh©n  corrections  ar©  applied,  th©  values  of  A  agree 
fairly  well  with  th©  atmospheric  pressure  values,  is  support  for 
at  least  the  approximate  validity  of  the  correction  procedure. 

The  variation  of  B  with  temperature  agrees  extremely  well 
with  the  atmospheric  density  case  for  wire  31.  It  is  probable  that 
this  close  an  agreement  is  fortuitous  in  view  of  the  experimental 
scatter-  The  uncorrected  values  of  B  show  a  slightly  different 
dependence  on  t+I  than  the  corrected  values  but  not  such  a  striking 
difference  as  for  the  quantity  A. 

Finally,  the  exponent  n  seems  to  vary  with  temperature  in  a 
way  somewhat  resembling  the  case  for  sma I ler  Knudsen  number;  again 
the  scatter  precludes  any  precise  comparison.  There  does,  however, 
appear  to  be  a  discernible  difference  in  the  absolute  value  of  n 
between  these  results  and  the  case  for  smaller  Knudsen  number.  The 
uncorrected  values  of  n  also  differ  considerably  from  the  corrected 
values,  the  former  giving  a  larger  dependence  on  wire  temperature. 

D.  Summary  of  Data 

For  ease  of  comparison,  a  summary  of  the  data  discussed  above 
is  given  in  the  table  below.  The  logarithmic  slopes  and  the  inter¬ 
cepts  at  T=0  are  given  for  each  case.  Included  are  the  results  from 
the  large  wire  at  0.45  atmospheres.  In  addition,  in  view  of  the  pre¬ 
vious  remarks  on  the  argon  results,  an  analysis  of  the  argon  data  was 
carried  out  assuming  a  temperature  Jump  correction  of  C  =  4  as  opposed 
to  the  value  C  =  2  taken  for  the  other  measurements.  For  all  cases, 
values  computed  without  corrections,  as  well  as  with,  are  also  given. 
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Table  2 


Pressure  Correction 

Wire  Gas  (atm)  Kn(To)  C^end  loss 


A  ‘  _B _ n_  A(0)  B(0)  n(0) 


30  Air 
Argon 

Air 

11,30  Air 

31  Air 


0.02  2,  E.L. 

0,0 

2,  E.L. 
0,0 

A ,  E.L. 


>2 

<0.02 

2,  E.L. 
0,0 

1 

0.02 

2,  E.L. 
0,0 

0.21 

0.04 

2,  E.L. 
0,0 

0.45 

0.02 

2,  E.L. 
0,0 

1 

0.01 

2,  E.L. 
0,0 

-0.04  0.15  -0.13  0.30  0.44  0.515 
-0.15  0.03  -0.20  0.32  0.45  0.49 
0.05  0.21  -0.10  0.23  0.50  0.45 
-0.21  0.14  -0.23  0.27  0.51  0.44 
-0.03  0.41  -0.09  0.27  0.48  0.49 
-0.27  0.33  -0.14  0.30  0.44  0.505 
-0.65  0.38  -0.18  0.31  0.47  0.50 
-0.04  0.14  -0.13  0.31  0.44  0.51 
-0.18  0.02  -0.22  0.33  0.45  0.49 
-0.07  0.36  -0.04  0.58  0.43  0.62 
-1.1  0.63  0.87  0.87  0.38  0.68 

-0.12  0.39  -0.02  0.58  0.40  0.60 
-0.60  0.55  -0.43  0.84  0.35  0.68 
-0.33  0.46  -0.26  0.64  0.43  0.59 
-0.29  0.18  -0.18  0.72  0.52  0.53 


The  following  additional  remarks  may  be  made  about  these  data: 

1)  The  differences  In  variation  of  A,  B,  n  with  t+I  between  corrected 


and  uncorrected  values  are  in  each  case  considerable  and  are  of  course 
most  noticeable  when  the  Knudsen  number  is  appreciable  (approximate 
values  of  Knudsen  number  based  on  free  stream  properties  are  given  in 


the  table  as  Kn(To))‘ 

2)  Only  two  significant  figures  are  given  for  each  quantity  since 
certainly  no  more  accuracy  can  be  attached  to  the  data  than  this, 
however,  in  view  of  interest  in  determining  the  exponent  n  for  the 

O 

case  T  -•  0,  the  rounding-off  was  done  to  three  figures  for  the  case 
of  air  at  atmospheric  pressure  (wire  30)  and  at  high  pressure.  It 
is  seen  that  the  exponent  may  be  taken  either  as  0.52  or  0.51  in 
the  first  case  and  0.50  or  0.51  in  the  second,  when  rounded  to  two 


f i gures . 

3)  Use  of  SC  =  4  for  the  argon  data  has  brought  the  results  for 
A(0),  B(0),  n(0)  closer  to  the  values  obtained  for  air,  although 
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tion  (ll)>  an  accomodation 


still  not  in  agreement.  From  equa 
coefficient  of  about  0.6  Is  implied  with  use  of  C  -  4  as 
opposed  to  f  =  0.95  for  C  =  2. 

E.  Heat  Loss  in  Helium 

A  striking  demonstration  of  the  temperature  jump  effect  is 
provided  by  several  heat  loss  measurements  carried  out  in  helium. 
Because  of  the  low  density  of  helium,  the  wire  Knudsen  number  for  the 
lx  10-^  in.  diameter  wires  Is  of  the  order  of  0.1  at  atmospheric 
pressure.  In  addition,  it  seems  to  be  well  established  that  the 
accomodation  coefficient  for  helium  on  material  such  as  platinum*  is 
considerably  smaller  than  for  air.  The  combination  of  these  factors 
results  in  a  large  value  of  the  temperature  Jump  correction  and  conse¬ 
quently,  without  the  correction,  a  marked  decrease  in  the  quantity  A  is 
to  be  expected.  Such  a  decrease  was  in  fact  found.  Figure  7  gives 
results  obtained  for  a  platinum  wire  of  lx  10"^  in.  diameter  operated 
at  a  fairly  low  value  of  T  in  air,  argon  and  helium.  The  lower  circles 
are  the  uncorrected  values  of  Nusselt  number  in  helium,  and  data  on 
the  upper  right  represent  the  Nusselt  numbers  recorded  for  air  and 
argon.  In  view  of  the  large  correction  required  for  helium  and  the 
lack  of  precise  knowledge  concerning  its  accomodation  coefficient, 
the  data  in  air  and  argon  were  used  to  compute  (on  the  basis  of 
equation  (II))  the  accomodation  coefficient  for  h6lium  and  therefore 
the  correction,  C-  Specifically,  at  a  single  Reynolds  number,  the 
corrected  Nusselt  number  for  air  or  argon  (which  was  essentially 
Identical  for  these  gases)  was  compared  with  the  uncorrected  value 

*0r  more  precisely,  for  helium  on  air  on  platinum.  The 
accomodation  coefficient  depends  strongly  on  the  surface 
condition,  e.g.  see  Kennard  (1938),  p.  320  et  seq. 


for  helium  and  the  correction  C  required  to  bring  the  helium  data 
into  coincidence  with  the  other  data  was  computed.  Then  this  value 
of  C  was  applied  to  the  remainder  of  the  helium  data  at  other 
Reynolds  numbers.  The  resulting  corrected  values  for  helium  are 
plotted  as  the  crosses  in  Figure  7.  For  these  data,  the  computed 

value  of  C  was  5.0  as  noted  on  the  figure. 

In  Figure  8,  results  from  a  different  wire  (of  the  same 
nominal  diameter)  at  higher  values  of  T  are  plotted.  Again,  the  lower 
data  are  for  helium  uncorrected  for  temperature  jump  and  the  same 
procedure  was  applied  in  computing  the  value  of  C  required  to  bring 
the  helium  data  into  colncid^ce  with  the  air  data  at  a  single  value 
of  Reynolds  number.  It  is  noted  that  a  value  of  C  of  11.7  was  ob¬ 
tained;  this  is  considerably  in  excess  of  that  found  at  the  lower 
value  of  T.  It  should  be  mentioned  that  these  values  of  C  were  com¬ 
puted  assuming  a  Knudsen  number  based  on  gas  properties  evaluated  at 
the  free  stream  temperature. 

The  different  values  for  C  found  at  different  wire  temperatures 
indicate  the  strong  dependence  of  accomodation  coefficient  on  wire 
temperature.  In  Figure  9  the  values  of  accomodation  coefficient  com¬ 
puted  (by  the  above-described  computation  process)  at  each  wire 
temperature  employed  are  plotted  against  film  temperature.  The  lower 
set  of  data  (unflagged)  corresponds  to  values  computed  by  taking  a 
Knudsen  number  based  on  free  stream  temperature,  KpfTo^^  and  the  upper 
set  (flagged)  by  basing  it  on  wire  temperature.  The  dotted  line  repre¬ 
sents  data  reported  by  Hartnett  (1961)  which  were  taken  by  Oliver  and 
Farber  for  helium  on  platinum  (they  were  obtained  by  a  conduction  heat 
loss  measurement  of  f ) .  The  agreement  in  trend  with  wire  temperature 
is  good  although  the  absolute  values  are  somewhat  different. 
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F.  Conclusions 


What  conclusions  may  be  drawn  concerning  the  heat  loss  from 
cylinders  at  low  Reynolds  numbers  from  the  above  information?  The 
specific  questions  in  point  are: 

1)  Can  the  heat  loss  law  be  accurately  represented  by  an 
equation  of  the  form 

Nu  =  A  +  BR  ? 

2)  If  so,  how  do  the  quantities  A,  B,  n  vary  with  temperature? 

Does  this  temperature  variation  depend  on  Knudsen  number? 

3)  Can  the  data  at  non-neg I i g i b  le  Knudsen  number  be  reduced 
to  the  case  of  Kn  -  0  by  the  simple  relation  (10)? 

With  regard  to  question  I ) ,  an  affirmative  answer  is  supported 
by  the  following  evidence:  for  a  given  wire  diameter  (wires  30  and  II) 

In  a  given  gas  (air),  nearly  identical  values  were  found  for  A,  B,  n 
at  T  -  0  whether  the  wire  was  operated  at  atmospheric  density  and  the 
velocity  alone  varied  or  whether  the  wire  was  operated  at  density 
levels  up  to  eight  times  atmospheric  and  the  density  and  velocity 
simultaneously  varied  In  addition,  for  a  wire  2.5  times  this  diameter, 
a  nearly  identical  value  for  B  at  t  -  0  was  found  over  a  range  of  den¬ 
sity  levels  from  atmospheric  to  about  0.2  atmospheres.  In  contradiction 
to  these  conclusions,  when  the  smaller  wire  was  operated  in  argon, 
somewhat  different  values  for  A,  B,  n  at  t  —  0  were  found;  further,  a 
significant  difference  from  the  small  wire  values  was  noted  in  the 
values  of  A  and  n  for  the  large  wire  over  all  density  levels.  There 
are  two  explanations  which  may  be  offered  to  account  for  these 
discrepancies,  over  and  above  the  experimental  scatter: 

First,  the  corrections  employed  may  not  be  strictly  applicable: 
this  premise  is  related  to  3)  above  and  will  therefore  be  discussed  later. 
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Second,  the  simple  heat  loss  taw  assumed  Is,  In  fact,  not 
precise  but  merely  a  rough  approximation  to  the  actual  law. 

It  appears  that  further  study  of  the  specific  discrepancies 
noted  is  required  to  decide  among  the  alternative  explanations  given. 

The  second  explanation  leads  directly  to  the  problem  of  what  other 
parameters  are  relevant.  The  discrepancies  in  the  large  wire  data 
might  be  attributed  to  an  aspect-ratio  effect  (the  possibility  of 
which  was  suggested  by  Cole  and  Roshko,  1954);  this  is  discussed  later. 
With  regard  to  the  argon  data,  the  only  possible  additional  parameters 
are  Pr  and  y  but  in  what  way  these  might  enter  in  the  heat  loss  law  is 
unknown;  the  difference  in  Pr  between  air  and  argon  is  only  a  matter 
of  a  few  percent  at  any  rate. 

Turning  to  question  2),  it  is  apparent  from  the  results  dis¬ 
cussed  that  A,  B,  n  all  vary  somewhat  with  wire  temperature.  In  addition, 
even  after  corrections  have  been  applied,  there  is  a  difference  between 
the  variation  with  temperature  at  non-neg I i gi b le  Knudsen  number  and  for 
Kp  -  0  at  least  in  the  quantities  A  and  B.  Thus,  while  the  temperature 
jump  correction  may  well  account  for  the  overall  decrease  in  Nusselt 
number,  it  does  not  apparently  take  into  account  the  effect  of  the  Jump 
on  the  temperature  distribution  about  the  wire  and  the  consequent  differ¬ 
ence  in  dependence  of  heat  transfer  on  wire  temperature  level  at 
different  Knudsen  numbers.  The  values  of  A  and  B  appear  to  decrease 
in  magnitude  as  the  Knudsen  number  increases  both  for  the  sma  I  I  and 
large  wires.  However,  n  seems  to  be  constant  for  the  sma  I  I  wires  and 
diminishes  somewhat  for  the  large  wire  with  increasing  Knudsen  number. 
There  is  insufficient  data  to  obtain  a  meaningful  plot  of  A,  B,  n  vs. 
Knudsen  number. 
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To  properly  answer  question  3)  requires  more  precise  informa¬ 
tion  on  the  accomodation  coefficient.  Since  the  Knudsen  numbers  in  a 
are  all  quite  small,  as  are  the  corresponding  corrections,  it  is  diff 
cult  to  make  accurate  conclusions  on  the  validity  of  the  temperature- 
jump  correction.  In  the  case  of  helium,  where  the  Knudsen  number  is 
considerably  larger,  the  accomodation  coefficient  is  not  adequately 
known.  Nonetheless,  the  general  effect  of  the  correction  on  the 
data  of  wire  31  as  shown  In  Figure  6,  as  well  as  the  fact  that  the 
application  of  the  correction  to  the  helium  data  appears  to  produce 
results  consistent  with  the  other  heat  loss  data  (and  gives  a  reason¬ 
able  value  for  the  accomodation  coefficient  and  Its  variation  with 
temperature),  all  indicate  the  general  validity  of  the  correction. 

The  improved  agreement  in  the  argon  data  when  a  sma I ler  accomodation 
coefficient  is  assumed  may  give  further  support  to  this  view,  pending 
the  determination  of  accurate  values  for  f. 

G.  Comparison  with  Results  of  Other  Investigators 

The  heat  loss  quantities  found  will  first  be  compared  with 
previous  results  for  the  case  t  -•  0.  Two  typical  theoretical  solu¬ 
tions  will  be  considered:  the  approximate  result  given  by  King 

15)  Nu  =  - ^ -  RPi"  <  0.08 

8 

jtn  f^pc  ~  0.67E 


16) 


RPR  >  0.08 


and  the  Oseen  solution  by  Cole  and  Roshko, 


17) 


Nu 


■r  -  0 


In  these  equations  E  is  Euler's  constant. 

Of  the  available  experimental  results,  the  most  precise 

appear  to  be  those  of  Co  1 1 is-W I i I iams.  They  give 
18)  Nu  =  0.24  +  0.56  0.02  <  R  <  44 


19) 


Nu 


2.09 

jfcn  1 1 .4/R 


0.02  <  R  <  0.5 


The  often  quoted  correlation  of  McAdams  (1954)  of  a  wide  variety  of 
data  from  various  investigators  is 


20) 


Nu  *  0.32  +  0.43  R®’^^ 


0.  I  <  R  <  1000 


The  result  obtained  in  the  present  work  for  wire  30  at 
atmospheric  pressure  in  air  (with  corrections  based  on  Kn(Tw) 
app lied)  is 

2n  Nu  =  0.30  +  0.44  R°'^^  0.4  <  R  <  4 


where  the  exponent  has  arbitrarily  been  rounded  off  from  0.515 
to  0.52.  In  argon, 

22)  Nu  =  0.23  +  0.50  0.4  <  R  <  4 

For  purposes  of  compar ison,. a  Prandt I  number  of  0.70  will  be 
taken  in  15)  to  17).  Each  of  the  equations  i5)  to  22)  is  plotted  in 
Figure  10  over  the  limits  specified.  With  the  exception  of  the 


-  27  - 


result  of  King  for  large  Reynolds  number,  there  appears  to  be  only 
small  differences  among  the  various  results.  As  pointed  out  by  Col  I  is 
and  Williams,  their  data  approach  the  theoretical  result  of  Roshko  and 
Cole  at  small  Reynolds  number.  At  higher  Reynolds  number,  the  Collis- 
Williams  data  appear  to  be  somewhat  above  the  results  of  the  present 
work  which  seems  to  agree  closely  with  McAdams'  correlation.  There 
is  only  a  very  small  difference  between  the  argon  and  air  curves  In 
the  present  work. 

It  has  been  mentioned  previously  that  Cole  and  Roshko  have 
suggested  a  possible  dependence  on  aspect  ratio  of  the  heat  loss  from 
finite  cylinders.  For  pure  conduction  from  an  ellipsoid  they  point 
out  that 

2 

23)  Nu  =  - —  (R  =  0> 

in2  R 

for  large  aspect  ratio.  Taking  /k  =  800  which  corresponds  to  the 
value  for  wire  30  then  from  23) 

Nu  =  0.27  (R  =  0) 

and  this  value  is  given  on  the  ordinate  axis  in  Figure  10.  It  should 
be  noted  that  this  value  is  very  close  to  the  value  of  A  obtained  for 
wire  30.  This  fact,  combined  with  the  finding  that  the  variation  of 
A  with  temperature  is  nearly  zero  (at  least  for  the  corrected  data 
obtained  at  atmospheric  pressure),  is  strong  indication  that  the 
quantity  A  in  the  assumed  heat  loss  law  represents  pure  conduction 
heat  loss  from  the  cylinder.  If  so,  then  it  may  well  be  aspect- 
ratio  dependent  in  a  way  quite  analogous  to  23).  The  somewhat 
different  variation  of  A  with  temperature  in  the  high  density  tests 
rema i ns  to  be  exp  I  a i ned ,  however . 
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about  300  hence  the  pure 


The  aspect-ratio  of  wire  31  was 
conduction  Nusselt  number  Is  0.31,  a  value  considerably  less  than  the 
value  of  0.64  obtained  for  wire  31  at  atmospheric  density.  It  does 
not  then  appear  possible  to  explain  the  discrepancy  In  A  between 
wires  31  and  30  on  the  basis  of  different  aspect-ratio. 

With  regard  to  the  heat  loss  law  for  finite  t,  a  wide 
divergence  in  its  dependence  on  t  appears  in  the  literature.  From 
a  consideration  of  many  factors,  the  results  of  Co  I  I is-Wi I  I  lams 
again  appear  to  be  the  most  precise.  Unlike  the  work  carried  out 
here,  however,  these  investigators  grouped  the  results  at  all  wire 
temperatures  together  to  determine  a  single  Reynolds  number  exponent, 
independent  of  wire  temperature.  For  this  reason.  It  is  difficult  to 
compare  their  results  with  the  present  results.  It  is  significant, 
however,  that  the  exponent  found  by  them  was  0.45,  a  value  less  than 
that  Which  has  been  reported  previously  in  the  literature.  The  work 
here  indicates  that,  in  general,  the  exponent  diminishes  with  wire 
temperature  from  a  value  near  0.5  at  T  -  0  to  lower  values  at  large 
temperatures.  Consequently,  grouping  together  results  taken  over  a 
range  of  wire  temperatures  would  naturally  tend  to  yield  an  exponent 
lower  than  0.5  as  the  best  data  fit  (this  fitting  was  done  by  "trial^ 
and  error"  In  the  work  of  Col  I  Is  and  W i I  I iams) .  It  is  bel ieved  then 
that,  for  the  assumed  heat  loss  law  given  by  equation  (I),  values  of 
the  exponent  n  lower  than  0.5  (as  found  by  Co  1  I i s-W i I  I i ams)  result 
from  the  effect  of  temperature  loading;  when  t  -  0,  the  exponent  is 
very  nearly  0.5,  which  is  the  conventionally  assumed  value. 
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VI.  VELOCITY  AND  TEMPERATURE  FLUCTUATION  SENSITIVITY 
A.  Pre I Imi nary  Remarks 

Once  the  steady  heat  loss  law  for  a  hot-wire  is  known,  the 
sensitivity  of  the  wire  to  fluctuations  in  the  stream  passing  over  it 
may  be  obtained  by  perturbation  of  the  mean  flow.  A  linearized  sensi¬ 
tivity  equation  thereby  results  for  the  (small)  fluctuation  quantities 
in  which  each  such  quantity  is  weighted  by  a  sensitivity  coefficient, 
these  coefficients  provide  a  measure  of  the  contribution  of  the  quan¬ 
tity  in  question  to  the  total  AC  signal  output  of  the  heated  wire.  A 
general  derivation  of  the  sensitivity  of  the  hot-wire  to  fluctuations 
in  the  stream  velocity  and  temperature  will  first  be  given  here;  the 
details  of  the  transient  response  including  lag  effects  will  be 
ignored*.  Subsequently,  simplified  forms  of  the  velocity  sensitivity 
for  special  cases  will  be  deduced  and  the  errors  introduced  by  their 
use  discussed.  Finally,  the  effects  of  end  loss  to  the  wire  supports 
will  be  considered.  Although  they  are  not  usually  important  in  hot¬ 
wire  anemometry,  for  purposes  of  estimation  the  effects  of  radiation 
and  a  non-linear  resistance-temperature  wire  behavior  will  be  in¬ 
cluded;  an  arbitrary  reference  temperature  upon  which  to  base  the 
fluid  properties  is  employed. 

B.  Sensitivity  Equation 

The  heat  loss  law  found  empirically  for  a  long  heated  cylinder 
placed  normal  to' an  incompressible  stream  at  low  Reynolds  number  and 
negligible  Knudsen  number  is  given  by  (I).  The  wire  resistance  is 

taken  as  ^  t 

Rw  =  Ro  r  I  +  Oo  (T^w”To)  +  Po  (Tw~To)  J 

*  For  a  discussion  of  this  response,  see  Corrsin  (1961). 
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consequent ly 


25) 


T  =  ( l“6) 

OLq 


to  order  where 


Using  a  film  temperature 

Tf  =  Tq  +  t  (Ty^~TQ) 

where  t  is  an  arbitrary  constant,  the  fluid  properties  evaluated 
at  this  temperature  are  taken  to  be 

26)  k  =  ko  (  I  +  2  thaw) 

27)  u  =  TJo  (  I  +  2  th|aw)  (1+2  thzaw) 


where 


h 


(dk/dT)^ 

- (1-6),  h, 


(d|x/dT)^ 

^^M'o 


1-6 

(1-6),  h2  =  - 


In  the  expression  (I),  a  functional  dependence 

A  =  A(t) 

B  =  B(t) 
n  =  n(T) 

is  assumed.  The  quantity  A  is  taken  to  Include  the  heat  loss  due  to 
radiation  so  that  a  separation 
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28) 


A(t)  =  Ar(t)  +  Af<T) 


is  effected  where  Ar  represents  heat  loss  due  to  radiation  and  Af  the 
value  of  A  in  the  absence  of  radiation.  For  the  radiation  heat  loss, 
the  Stefan-Boltzmann  law  is  applied  and  the  wire  assumed  to  radiate  to 
a  large  enclosure  which  is  at  the  ambient  stream  temperature. 
Consequent iy. 


29) 


Ar  = 


J  T  4  T  4 
doe  “Tq 


k  Tw-Tc 


where  a  is  the  Stefan-Boltzmann  constant  and  e  the  wire  emissivity. 

The  fluctuation  quantities  of  the  field  relative  to  their  mean 
value  are  denoted  by  AU/U  and  ATq/To;  these  lead  to  a  fluctuating  output 
voltage  of  the  DC  heated  wire,  AE/E.  This  voltage  may  be  expressed  as 


30) 


^  ^  ^l_JRw  j 

E  Rw  V 


where 


31) 


V  =  -|2 


ARw 

Al 


and  is  determined  from  the  characteristics  of  the  circuit  into  which 
the  wire  is  placed.  In  the  usual  constant-current  hot-wire  circuit 
where  the  wire  forms  part  of  a  bridge  which  is  placed  in  series  with 
a  large  resistor  and  heated  from  a  battery, 

V  =  Eb  . 


the  battery  voltage. 
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To  obtain  the  fluctuation  sensitivity  equation,  the  expressions 
(24)  to  (29)  are  inserted  into  (I)  and  small  perturbations  of  U,  it 

are  assumed.  Linearizing  with  respect  to  these  quantities  and  using 
(30)  then  gives  an  expression  of  the  form 


32) 


AU  ^  ATo 

U  2  To 


S3 


E 


where  S j ,  S2  are  then  the  velocity  and  temperature  fluctuation 
sen si t I vit les. 

Denote  by  a  the  logarithmic  derivative  of  a  quantity  with 
respect  to  t+I ;  thus  e.g.. 


d/nA  , 

djtn  (t+I) 


also  call 


The  velocity  sensitivity  is  then  given  by 

A 

33)  S|  =  ”aw*^^^”Nu^ 

and  for  the  temperature  sensitivity. 


34) 


^2  2h2 


where 


■^o 


2thav»  p  A _ 

l+2thaw^  Nu  ’’o 


+  €0 


]  . 


In  this  expression. 
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wherel n 


\  Wn  ^  (  I+t)"^  -  (  H4t)  t 

[  I  .s,  - ^ - ] 


I  n  l+2thaw  ^^2  ^ 


2tt  I  +h2®w  ^ 


P  =  I  +  So  a„(l42thaw) 


h2  l+2thav<  p  <j-+l  ^ 

“b  =th  l+2h2a^  '  "17“  ^  ^ 


‘+^2aw  ill  .  W  r  ,  _  T+l.  .  -| 

=  1  +  2tT  1 42t  h I ^2  n  T+l  hr 


Afko(l-6) 


i  -  6  h2 
”  uto++ 


w  =  B  +  n  jfn  -  I) 


The  quantity  S3  is  given  by 


I  +  2a^-^  +  (1+avy)  f| 
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where 


2thaw  A  ,  ,  A  - 

f  =  -  —  r ,  +  e ,  (  I  -  7377) 

1 42thaw  Nu  I  •  Nu 


Pi  =  ~  (  I  +0)1  ) 

I  n  ' 


I  +S2 


I-(I+t)^  (I-3t) 


e I  =  (1-8)  €0 


(Di  =  Ofe  11“  .  6) 

'  ,  Ufe 


S2  =  <  l“<)  1-0),  S| 


C.  Simplifications  for  Velocity  Sensitivity 

It  is  clear  from  these  results  that  a  considerable  amount  of 
Information  is  required,  to  compute  the  sensitivities  in  general.  An 
alternate  method  of  desfr>ermining  the  fluctuation  quantities  in  the 
field  Is  the  so-called  "mode-diagram"  method  of  Morkovin  (1956).  In 
this  method)  a  form  of  the  heat  loss  law  is  assumed  simi  lar  to  (  I)  but 
with  n  =  i.  Then  instead  of  using  analytical  expressions  as  given 
here)  a  graphical  solution  technique  is  developed.  The  data  required 
for  this  solution  are  a  number  of  values  of  the  heat  loss  and  AC  wire 
output  signals  at  different  wire  temperatures;  having  this  information 
then  permits  a  I  I  of  the  fluctuation  quantities  to  be  computed.  The 
simplicity  of  the  mode— diagram  method  rests  upon  the  assumption  of  a 
fixed  exponent  of  i  for  the  Reynolds  number  in  the  heat  loss  law. 
However)  in  the  work  here)  the  exponent  was  found  to  be  a  fund  ion 

O 

of  wire  temperature.  Another  objection  to  the  graphical  solution 
method  is  that  a  large  number  of  wire  output  signals  must  be 
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obtained;  since  for  each  of  these,  a  separate  square-wave  compensa¬ 
tion  adjustment  must  be  made  each  of  which  entails  several  bridge 
balancing  operations,  the  method  is  time-consuming  and  is  not  practical 
when  running  time  is  limited.  Further,  the  method  precludes  obtaining 
Information  on  the  wire  heat  loss  under  various  flow  conditions  which 
Is  useful  in  checking  that  the  wire  is  behaving  properly. 

With  use  of  the  constant-cixrent  fluctuation  measurement  tech¬ 
nique,  a  direct  calibration  procedure  in  which  the  current  is  held 
fixed  and  the  wire  temperature  vs.  velocity  characteristic  is  obtained 
might  also  be  employed  and  is.  In  fact,  the  most  logical  method. 
However,  the  form  (1)  indicates  that  considerable  difficulty  would 
be  encountered  in  empirically  fitting  data  obtained  in  this  fashion. 

When  velocity  fluctuations  are  of  prime  interest,  a  direct 
analytical  computation  employing  the  general  information  previously 
obtained  for  the  heat  loss  law  may  be  applied  when  the  wire  Is  oper¬ 
ated  so  as  to  cause  it  to  be  principally  sensitive  to  velocity 
fluctuations.  The  following  simplifications  may  then  usually  be 

made: 

^  •  X  ■  _  n 

1)  No  free  stream  temperature  variations,  - —  - 

*o 

2)  Effects  of  radiation  negligible;  this  amounts  to  taking 
zero  emissivity  hence  Sq  =  S|  =  52  =  0- 

3)  Linear  resistance-temperature  behavior,  6=0. 

With  these  assumptions,  and  taking  f  =  i,  the  quantity  f|  in 


S3  becomes 


hay,  A 
I  +hay,  Nu 


l+«)|  ^  ®  1 


where  now 
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21^  14-ha,^  * 

h  I  +t>23w 


«l 


l-n 


h2  l+ha^  g 

h  I  +^'2®w 


h  I  I +^2^w  2w  I +h23w 

9=1+ - - - , 

h2  l+h|aw  l+2h2aw 

-  -  ,  A  r  Nu  -  II 

w=B+nXng  J 

and  S,  remains  the  same.  The  sensitivity  equation  then  is 

AU  _  S3  AE 
U  S,  E 

and  to  compute  the  velocity  fluctuation  from  the  wire  output 
siqnal  requires  a  knowledge  ofJ 

a)  the  fluid  properties  ko,  (dk/dT)m>  1*0^  (dn/dT)^, 

b)  the  operating  Nusselt  number,  Nu,  and  overheat  ratio,  a^ 

c)  the  heat  loss 'bonstants"  A,  B,  n  at  the  wire  temperature 
employed  and  their  variation  with  temperature,  A,  B,  n 

d)  the  coefficient,  and  stream  temperature,  T© 

e)  the  voltage  across  the  wire,  IRw,  and  the  supply  voltage,  V. 

The  degree  of  approximation  made  in  assumptions  1)  to  3)  will 

be  estimated. 

1)  Assumption  of  ATq/Tq  =  0: 

To  estimate  the  effect  of  temperature  fluctuations,  the 
assumptions  2)  and  3)  are  made  (for  simplicity)  so  that  S2  becomes 

S2  =  ^  (i-fo)  ; 
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hence 


The  ratio  S2/S |  provides  a  measure  of  the  relative  contribution  to 
the  total  voltage  signal  of  the  temperature  and  velocity  fluctuations. 

It  is  clear  that  within  the  range  of  allowable  wire  temperatures,  the 
wire  is  principally  sensitive  to  temperature  fluctuations.  Inasmuch  as 
It  is  not  possible  to  reduce  the  relative  temperature  sensitivity  to 
zero,  the  contribution  of  small  stream  temperature  fluctuations  to  the 
wire  output  signal  must  be  assessed.  For  an  Isotropic  field,  the  corre' 
lation  of  the  (scalar)  temperature  with  the  velocity  field  is  zero;  the 
less  restrictive  assumption  of  axial  symmetry,  however,  allows  a  non¬ 
zero  correlation  but  since  no  rational  choice  is  available  for  the 
correlation  in  this  or  in  a  more  general  field,  the  field  will  be 
assumed  to  be  isotropic  for  the  purpose  of  the  present  estimates. 

Then,  the  ratio  of  actual  rms  velocity  fluctuation  to  that  computed 
by  assuming  no  temperature  fluctuation  may  be  shown  to  be,  for  small 
temperature  fluctuations, 

U’act  '  ^ 

'"2  \S||  I  U'/U  1 

To  obtain  an  estimate  of  orders  of  magnitude,  assume 

36) 

and  take  the  typical  values  A/Nu  =  i,  n  =  i;  then 

S2  _  £ 

s7  ~  T  2+t 
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Examples  of  the  errors  introduced  with  sina  I  I  and  large  temperature 
fluctuations  (about  0.01  and  0.  I^C.  in  absolute  value)  when  the  wire 
is  operated  at  small  and  large  values  of  t  for  a  range  of  relative 
rms  velocity  fluctuations  are  given  in  the  following  table: 

Table  3 

T  Tj>/To  i%)  UVU  Error  (f) 


0.003 

0.01 

>100 

0. 1 

70 

1  .0 

0.7 

0.03 

0.01 

>100 

0.  1 

50 

1  .0 

0.5 

In  conventional  flow  streams  produced  in  the  laboratory,  a  level 
of  temperature  fluctuations  in  excess  of  0.0l°C.  may  easily  occur  unless 
very  special  precautions  are  taken.  From  the  above  table.  It  is  then 
clear  that,  unless  high  wire  temperatures  are  employed,  the  measurement 
of  small  velocity  fluctuations  are  at  best  only  approximate  and  may  be 
grossly  in  error.  Even  when  the  wire  is  operated  at  600°K.,  the  error 
in  measuring  velocity  fluctuations  of  the  order  of  0.1^  is  about  50f 
when  temperature  fluctuations  of  only  0.0l°C.  occur.  Consequently,  the 
accuracy  of  measurements  reported  in  the  literature  of  such  small  vel¬ 
ocity  fluctuation  levels  is  questionable  if  some  indication  that  the 
temperature  fluctuation  level  is  negligible  is  not  given.  A  simple 
way  to  provide  such  an  indication  experimentally  is  to  operate  the  wire 
at  more  than  one  (high)  temperature  level  and  compute  the  turbulence 
level  from  each  measurement,  assuming  no  temperature  fluctuation.  If 
the  resulting  values  agree  closely,  then  the  measurements  are  probably 
sat i sfactory . 

Since  it  is  necessary  to  operate  at  high  wire  temperature  to 
minimize  the  contribution  of  possible  stream  temperature  fluctuations 
to  the  wire  output  signal,  a  knowledge  of  the  wire  heat  loss  law  for 
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such  elevated  temperatures  Is  mandatory.  That  Is,  terms  due  to  the 
quantities  A,  B,  n  become  more  important  as  t  Increases.  As  an 
example,  consider  the  contribution  of  A  alone.  Making  the  assumptions 
of  (36),  and  again  taking  n  =  i,  Nu/A  =  i.  It  is  a  simple  matter  to 
show  that  velocity  fluctuations  computed  by  including  A  relative  to 
the  values  computed  assuming  A  is  zero  are,  for  small  T,  equal  to 


and  consequently  the  error  committed  in  ignoring  A  increases  directly 
with  T.  Since  A  may  typically  have  a  magnitude  of  0.2  then  at  t  =  i 
the  error  is  of  the  order  of  and  will  be  larger  for  higher  values 
of  T.  Measurements  reported  in  the  literature  taken  for  high  values 
of  T  may,  therefore,  also  be  considerably  in  error  If  the  temperature 
dependence  of  the  quantities  A,  B,  n  is  not  determined.  It  should  be 
noted  further  that  since  the  computation  of  the  turbulence  energy 
level  Involves  the  square  of  the  velocity  fluctuation,  the  errors 
committed  are  then  approximately  doubled. 


2)  Assumption  of  negligible  radiation 

The  effect  of  wire  radiation  is  determined  by  the  magnitude 

A 

of  the  quantities  Sq,  S|,  S2.  Taking  t  =  i,  A  =  0,  6=0,  and  using 
(36)  gives 

S2  =  s  I  =  2So 


where 


Sr>  = 


^f’^o 
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Using  »  lib.r»l  vsin.  of  0.3  for  th.  Wilssivlty  =nd  typicol  voliios 
of  0.25  for  Af,  300OK.  for  T,,,  gives  'O'-  « "-es  of  2.5*10-*  inches 
diometer  in  oir  o  voice  of  s„  of  oboct  SxiO-*.  Thus,  the  contribu¬ 
tion  of  the  rodlotlon  terms  is  less  then  \f  for  the  moximum  values 
of  T  ehlch  may  be  employed  and  these  terms  may  then  safely  be 
omitted  even  for  hires  of  relatively  large  diometer. 

3)  Assumption  of  6  =  0 

The  difficulty  In  precisely  measuring  the  second  coefficient 
makes  the  estimation  of  6  uncertain.  Spangenberg  (1955)  has  given 
several  values  of  go  90Pt-l0Rh  wires  <hich  differ  considerably. 

The  maximum  value  of  Po/Oo^ 

minimum  value  Is  zero.  A  simple  average  of  these,  say  -0.05  might 
be  a  fair  estimate  in  general.  Taking  this  value,  then  for  a^  =  0.5, 

6  becomes  0.025  hence  an  error  of  about  2.5^  may  be  Incurred  by  taking 
6=0.  it  should  be  pointed  out,  however,  that  if  6  =  0  Is  also 
assumed  in  the  heat  loss  data  reduction  then  the  use  of  emp i r ica 1 ly 
determined  quantities  A,  B,  n  will  partially  counteract  errors  in 
the  sensitivity  expression  which  results  from  taking  6  =  0.  The 
ultimate  test  of  course  is  whether  the  same  fluctuations  are  computed 
when  the  wire  is  operated  at  different  temperatures. 

D.  Effect  of  Wire  End  Losses 
I,  On  computing  fluctuation  level. 

A  consideration  of  the  effect  of  heat  loss  to  the  wire  supports 

on  the  measurement  of  fluctuation  quantities  has  been  given  by  Betchov 
(1952).  The  results  there  are  given  in  terms  of  the  (small)  parameter 

S  (see  equation  (9)). 
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The  theory  of  Betchov  was  carried  out  assuming  a  heat  ioss 
iaw  of  the  form 

Nu  =  A(t)  +  B 

where  Nu  and  R  are  based  on  stream  temperature^  A  Is  a  function  of 
wire  temperature  and  B  a  constant*  it  is  therefore  an  extension  of 
earlier  theories  in  which  A  Is  taken  as  a  constant,  although  it  does 
not  account  for  possible  variations  in  B  and  n  with  wire  temperature. 
It  may  be  shown,  however,  that  to  order  S  the  effect  of  end  ioss  and 
of  variable  A  are  separable  and  that  consequently  end  ioss  corrections 
may  be  applied  independent  of  the  effect  of  variable  A.  It  is 
therefore  reasonable  to  suppose  that  thesame  situation  will  hold  with 
n  and  B  also  variable,  especially  since  A  Is  the  principal  variable 
of  the  three. 

Expanding  the  result  of  Betchov  in  powers  of  S  gives  for  the 
end- loss  correction  to  the  fluctuation  measurement, 

l-a^ 

K^l  =  I  +  i  ^  5  > 

ignoring  terms  of  order  and  smaller.  For  S  of  about  0.05,  then 
at  a^  =  0.5  the  correction  is  only  of  the  order  of  ]%,  This  estimate 
then  indicates  that  the  correction  may  usual iy  be  safely  ignored.  it 
should  be  noted  that  it  is  assumed  in  the  above  that  the  quantities 
A,  B,  n  which  are  used  are  those  for  the  actual  wire  heat  loss,  and 
are  not  values  corrected  for  end  loss  or  temperature  Jump. 

2.  On  time  constant 

The  question  of  the  effect  of  end- loss  on  the  time  constant 
of  the  wire  was  also  considered  by  Betchov.  it  was  found  that  the 
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wire  frequency  characteristic  is  approximatei y  the  same  with  or 
without  end- loss  and  the  required  correction  to  the  time  constant 
was  computed  (however,  see  Corrsin,  1961,  for  a  discussion  of  the 
applicability  of  this  computation).  Again,  it  may  be  shown  from 
these  results  that  the  effects  of  variable  A  and  end- loss  can  be 
separated  with  the  result  that  the  relation  between  the  time  constant 
with  and  without  end- loss  is,  to  order  S, 

M  _ _  S 

—7=1-  J  l+Qw  2 

which  indicates  an  expected  deviation  o.  only  a  few  percent.  Since 
the  time  constant  is  usually  determined  experimentally  by  the  square- 
wave  method,  it  is  then  necessary  to  know  the  effect  of  end-loss  on 
the  current-fluctuation  time  constant.  To  order  S,  the  current  time 
constant  is  identical  with  or  wi thout  end- loss,  consequently,  the 
value  determined  for  the  velocity  will  be  in  error  by  the  (small) 
amount  given  by  the  above  expression  for  the  case  considered  by 

Betchov . 

E.  Variation  of  Time  Constant  with  Wire  Temperature 

Without  the  end-ioss  effect,  the  result  for  the  time  conslant 

given  by  Betchov  may  be  expanded  to  give 

M  ^  24a  y 

38)  —  =1-  0  4- - H  GA 

Mo  l+aw 

ignoring  terms  of  order  g3  and  smaller.  In  this  expression,  M  is 
the  time  constant  for  A  =  A(t),  Mo  the  time  constant  for  constant  A, 
and  the  parameter  G,  which  is  a  measure  of  the  effect  of  variable  A, 
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a  Iso 


G  =  h2aw  ( l+Ow)  » 

Cy^Sy^d^  /  |4Qm\ 

=  4k  (  Nu  I 

Equation  (38)  is  plotted  in  Figure  II  for  two  values  of  ay*.  The 
exper imentij^  po I nts  were  taken  at  random  from  time  constant  values 
obtained  by  the  square  wave  method  in  the  course  of  other  work  per¬ 
formed  here  which  is  to  be  reported  on  later.  To  eliminate  the 
effect  of  uncertainties  in  the  wire  diameter  and  other  physical 
properties  of  the  wire,  a  normalizing  procedure  was  used  in  p iott i ng 
these  data.  For  a  given  wire,  the  value  of  M/K^  recorded  for  the 
lowest  value  of  G  was  divided  into  the  rest  of  the  experimental  values 
of  M/Mq.  Then,  since  this  value  of  G  was  not  actually  zero,  a  correc¬ 
tion  was  applied  by  finding  the  theoretical  value  of  M/Mq  sf  this  value 
of  G  and  correcting  the  remainder  of  the  ratios  by  multiplying  by  the 
theoretical  value  there.  Although  there  is  some  scatter  in (the  data, 
the  general  shape  of  the  theoretical  prediction  appears  to  be 
verified.  It  should  be  noted  that  since  only  terms  up  to  order  G^ 
were  included,  the  theoretical  curve  plotted  may  be  in  error  by 
several  percent  for  G  =  0.3  or  greater. 
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12.  Abstract:  An  experimental  study  of  the  forced  convection  heat 
loss  from  long  thin  cylinders  placed  normal  to  an 
incompressible  flow  has  been  carried  out.  Tests  were  made  In 
three  different  gases  (air,  argon,  helium)  at  atmospheric  pressure 
and  in  air  over  a  pressure  range  0.2  to  8  atm.  Particular  attention 
was  given  to  determining  the  effect  of  cylinder  temperature  on  the 
heat  loss.  It  was  found  that  the  "constants"  employed  in  the  usual 
empirical  representation  of  this  heat  loss  vary  with  cyl inder 
temperature  in  a  way  dependent  on  Knudsen  number,  results  in 
helium  giving  a  marked  indication  of  this.  A  velocity  fluctuation 
sensitivity  equation  was  derived  which  includes  the  effects  of 
this  temperature  variation. 


Figure  I.  Hot  wire  probe  tip 


Figure  2.  Actual  wires  used.  Optical  comparator  photograph, 
upper  wire  is  1x10“^  inches  diameter  (wire  30), 
lower  wire  is  2.5x10“'*  Inches  diameter  (wire  31) 
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FIG.  3  A  -  USUAL  DIMENSIONAL  HOT  WIRE  CALIBRATION 
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FIG.  3  B  -  DIMENSIONLESS  HEAT  LOSS  FOR  ABOVE  DATA. 
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WIRE  31  . 


IV  Kl- 


NUSSELT  NUMBER  ,  Nu 


10 


5 


1.0 


0.5 


0.1 


0.1 


WIRE  5  (Pt) 

r«0.23 


TEMPERATURE 
JUMP  CORRECTION 
CONSTANT^  C 

2.0 

O 

5.0 

2.0 


GAS  SYMBOL 

ARGON  ^ 
HELIUM  6 
HELIUM  + 

AIR  V 


jF 


+ 

-O- 


+ 

o 


+ 

o 


HELIUM  DATA  UNCORRECTED 
FOR  TEMPERATURE  JUMP. 


0.5 

REYNOLDS 


1.0 

NUMBER,  R 


FIG.  7  -  Nu  vs.  R  FOR  AIR,  ARGON,  HELIUM  -  WIRE  5 
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FIG.  8  -  Nu  VS.  R  FOR  AIR  AND  HELIUM  -  WIRE  16 
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FIG.  9  -  COMPUTED  ACCOMODATION  COEFFICIENT 

FOR  He  ON  Pt. 
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FIG.  1 1  -COMPARISON  OF  MEASURED  TIME 
CONSTANTS  WITH  THEORY  OF  BETCHOV. 
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